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We study both experimentally and theoretically the losses 
induced by parametric excitation in far-ofT-resonance optical 
lattices. The atoms confined in a ID sinusoidal lattice present 
an excitation spectrum and dynamics substantially diflterent 
from those expected for a harmonic potential. We develop a 
model based on the actual atomic Hamiltonian in the lattice 
and we introduce semiempirically a broadening of the width 
of lattice energy bands which can physically arise from in- 
homogeneities and fluctuations of the lattice, and also from 
atomic collisions. The position and strength of the parametric 
resonances and the evolution of the number of trapped atoms 
are satisfactorily described by our model. 
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I. INTRODUCTION 

The phenomenon of parametric excitation of the mo- 
tion of cold trapped atoms has recently been the sub- 
ject of several theoretical and experimental investigations 
iQ-H . The excitation caused by resonant amplitude noise 
has been proposed as one of the major sources of heating 
in far-off-resonance optical traps (FORTs), where the 
heating due to spontaneous scattering forces is strongly 
reduced In particular, the effect of resonant excita- 
tion is expected to be particularly important in optical 
lattices, which usually provide a very strong confinement 
to the atoms, resulting in a large vibrational frequency 
and in a correspondingly large transfer of energy from 
the noise field to the atoms |l| . 

Nevertheless, parametric excitation is not only a source 
of heating, but it also represents a very useful tool to 
characterize the spring constant of a FORT or in general 
of a trap for cold particles, and to study the dynamics of 
the trapped gas. Indeed, the trap frequencies can be mea- 
sured by intentionally exciting the trap vibrational modes 
with a small modulation of the amplitude of the trap- 
ping potential, which results in heating or losses [^|J|] 
for the trapped atoms when the modulation frequency 
is tuned to twice the oscillation frequency. This proce- 
dure usually yields frequencies that satisfactorily agree 
with calculated values, and are indeed expected to be 
accurate for the atoms at the bottom of the trapping po- 
tential. From the measured trap frequencies is then pos- 
sible to estimate quantities such as the trap depth and 
the number and phase space densities of trapped atoms. 



We note that this kind of measurement is particularly 
important in optical lattices, since the spatial resolution 
of standard imaging techniques is usually not enough to 
estimate the atomic density from a measurement of the 
volume of a single lattice site. 

Recently, ID lattices have proved to be the proper envi- 
ronment to study coUisional processes in large and dense 
samples of cold atoms, using a trapping potential inde- 
pendent for the magnetic state of the atoms. In this 
systems, the parametric excitation of the energetic vi- 
brational mode along the lattice provides an efficient way 
to investigate the cross-dimensional rethermalization dy- 
namics mediated by elastic collisions 

Most theoretical studies of parametric excitation rely 
on a classical [^j or quantum harmonic approximation 
of the confining potential. Under certain circumstances 
these expressions show quite good agreement with exper- 
imental results 1^. However, general features of the opti- 
cal lattice could be lost in these approaches. For example, 
a sinusoidal potential exhibits an energy band structure 
and a spread of transition energies, while harmonic oscil- 
lators have just a discrete equidistant spectrum. Thus, 
we might expect that the excitation process may happen 
at several frequencies, and with a non-negligible band- 
width. Such anharmonic effects can be important when- 
ever the atoms are occupying a relatively large fraction 
of the lattice energy levels. The purpose of this paper is 
to give a simple description of parametric excitation in 
a sinusoidal ID lattice. In the next Section, we briefly 
discuss general features of the stationary states on such 
a lattice. Then, we summarize the harmonic description 
given in Ref. j^] and extend it to the anharmonic case. 
By a numerical evaluation of transition rates, we make 
a temporal description of parametric excitation which is 
compared with experimental results. We discuss about 
the relevance of broadening of the spectral lines to un- 
derstand the excitation process in this kind of systems. 
Some conclusions are given in the last Section. 



II. STATIONARY STATES OF A SINUSOIDAL 
OPTICAL LATTICE. 



The Hamiltonian for an atom in a red detuned FORT 



IS 



(2.1) 
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with 



Kff(x) = --a|f(x)|2, 



(2.2) 



where a is the effective atomic polarizations and £{x) is 
the radiation field amplitude. For the axial motion in a 
sinusoidal ID lattice we can take 
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The corresponding stationary Schrodinger equation 
can be written in canonical Mathieu's form 
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It is well known that there exists countably infinite sets 
of characteristic values {a^} and {br} which respectively 
yield even and odd periodic solutions of Mathieu equa- 
tion. These values also separate regions of stability. In 
particular, for g > the band structure of the sinusoidal 
lattice corresponds to energy eigenvalues between and 
hr+i [0. The unstable regions are between hr and ar- 
For g >> 1, there is an analytical expression for the band 

width inli: 
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The quantities defined above can be expressed in terms 
of a frequency ujq defined in the harmonic approximation 
of the potential 
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thus obtaining 



V? 



T^) ■ (2-10) 



Thus, the width of the r-band can be estimated using 
Eq. ( p.8D whenever the condition {Vq/Tiujq)'^ » 1 is sat- 
isfied. In the experiment we shall be working with a ID 
optical lattice having Vq ~ 10.5?ia;o- While the lowest 
band r = has a negligible width ~ lO^^^TiuJo, the band 
widths for highest lying levels r = 10, 11, 12, and 13 would 
respectively be 0.0065 ,0.1036, 1.52, and 20.56 in units of 

hujQ. 



In order to determine the energy spectrum, a varia- 
tional calculation can be performed. We considered a 
harmonic oscillator basis set centered in a given site of 
the lattice, and with frequency ojq. The d iago nalization 
of the Hamiltonian matrix associated to (2.4) using 40 
basis functions gives the eigenvalues £'„ < Vq shown in 
Table I for Vq — W.^Hluq. According to the results of last 
paragraph, the eigenvalues 12 and 13 belong to the same 
band while the band width for lower levels is smaller than 



III. PARAMETRIC EXCITATION 

As already mentioned, parametric excitation of the 
trapped atoms consists in applying a small modulation 
to the intensity of the trapping light, 
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V,s[l + eit)]. 



(3.1) 



Within first order perturbation theory, this additional 
field induces transitions between the stationary states n 
and m with an averaged rate 



where 
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T(m, n) = {m\Vcs\n) 

= EnSr, 
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(3.3) 



is the matrix element of the space part of the perturba- 
tion and 



dr cos UJT{e{t)e{t + t)) 



(3.4) 



is the one-sided power spectrum of the two-time correla- 
tion function associated to the excitation field amplitude. 

If the confining potential is approximated by a har- 
monic well, the transition rates different from zero are 
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^^S'(0)(2n-|- 1) 



Rn±2^n = -^S{2ujQ){n + 1 ± 1)(71 ± 1) 
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(3.5) 
(3.6) 



The latter equation was used in to obtain a simple 
expression for the heating rate. 



{E) = -ujIS{2uq){E) , 



(3.7) 



showing its exponential character. The dependence on 
2ljq is characteristic of the parametric nature of the exci- 
tation process. The fact that fi is not present is consistent 
with the apphcability of Eq. (3.7) in the classical regime. 
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Classically, parametric harmonic oscillators exhibit 
resonances not just at 2wo but also at ^io^jn with n any 
natural number In fact, the resonances corresponding 
to n=2, i.e. at an excitation frequency cjq, have been ob- 
served in optical lattices |^,^ . A quantum description of 
parametric harmonic excitation also predicts resonances 
at the same frequencies via n-th order perturbation the- 
ory [p^ . In particular, the presence of the resonance at 
Wo can be justified with the following argument. Accord- 
ing to the standard procedure, the second order correc- 
tion to the transition amplitude between states |n) and 
|m) is given by 

i?(^L„ = {n\U'^^\t^,t)\m) = ^ (^)V(n,fc)r(fc,m) 

k 
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(3.9) 



These states are ordered according to their energy: iJ„ < 
i?„+i as exemplified in Table I. Within this scheme one 
obtains a very simple expression for T(n, m) 



T(n, m) = Enbn 



^n,^MP''\3). (3.10) 



I dt'e^^-^^'eit') [ dt" e'"^""^'" e{t") (3.8) 

Jta Jtn 



I to 

with U'^^\to,t) the second order correction to the evo- 
lution operator U. Therefore, the transition may be de- 
scribed as a two step procedure |m) <— \k) <— For 
harmonic parametric excitation the matrix element of the 
space part of the perturbation differs from zero just for 
transitions \n) <— \n) and |n±2) <— \n). Consider a tran- 
sition in Eq. ( |3.8D involving a "first" step in which the 
state does not change |n) ^ \n) and a "second" step for 
which |?i ± 2) <— \n). Then resonance phenomena occur 
when the total energy of the two excitations, 2hn, coin- 
cides with that of the second step transition, i. e. for an 
excitation frequency Q — lUq. 

These ideas can be directly extended to anharmonic 
potentials: the corresponding transition probability rates 
R{n, m) would be determined by the transition matrix 
T(n,m), by the transition frequencies Wnm and by the 
time dependence of the excitation e{t). In general, an- 
harmonic transition matrix elements T{n, m) will be dif- 
ferent from zero for a wider set of pairs (n^m). Besides, 
the transition energies will not be unique so that the exci- 
tation process is not determined by the excitation power 
spectrum at a single given frequency 2luq and its sub- 
harmonics iLOojn. As an example the transition energies 
for the specific potential considered in this work are re- 
ported in Ta ble I. Therefore, within the model Hamilto- 
nian of Eq. (3.1), resonance effects can occur for several 
frequencies that may alter the shape of the population 
distribution within the trap. However, in general these 
resonant excitations will not be associated with the es- 
cape of trapped atoms. 

Here we are interested in a ID lattice; the direct ex- 
tension of the formalism mentioned above requires the 
evaluation of the matrix elements T(n, m) among the 
different Mathieu states that conform a band. This in- 
volves integrals which, to our knowledge, lack an ana- 
lytical expression and require numerical evaluation. As 
an alternative, we consider functions which variationally 
approximate the Mathieu functions. They are the eigen- 
states of the Hamiltonian (2.4) in a harmonic basis set of 
frequency uji^: 



It is recognized that any discrete basis set approxima- 
tion to a system with a band spectrum will lack fea- 
tures of the original problem which have to be care- 
fully analyzed. Anyway, alternatives to a discrete ba- 
sis approach may be cumbersome and not necessarily 
yield a better approach to understand general proper- 
ties of experimental data. While the discrete basis ap- 
proach is exact for transitions between the lowest levels, 
which have a negligible width, eigenstates belonging to 
a band of measurably width should be treated with spe- 
cial care. Thus, we shall assume that matrix elements 
T(i/, /i) involving states with energies Ei, and f^, so that 

£;„ - (Er, - Er,-l)l2 <E^<Er, + (S„+l - Er,)/2 with 

an analogous expressions for are well approximated 
by T{n, m). 

Within this scheme the equations which describe the 
probability P{n) of finding an atom in level n, given the 
transition rates Rm^n are 



P(n)=^i?(^L„(P(m)-P(n)) 



(3.11) 



in the first order perturbation theory scheme, and the 
finite difference equations 

Pn{t) = Pnito) + 

Y,R[ltn{Pm{to) - PnitoM - to) + 

m 

^i?(^L„(P„(to) - Pnitomt - tof , (3.12) 

m 

valid up to second order time-dependent perturbation 
theory whenever t ^ to. Both sets of equations are sub- 
jected to the condition 



^P(n) = l. 



(3.13) 



(r) 

tion of Rn<—m also requires the specification of the spec- 
tral density S(lu). In the problem under consideration, 
the discrete labels m, n are used to calculate interband 
transitions which are actually spectrally broad. This 
broadening might arise not only from the band structure 
of t he e nergy spectra associated with the Hamiltonian 
Eq. ( ^.4[ ), but also from other sources, which we will dis- 
cuss below. Broad spectral lines can be introduced in our 



Now, according to Eqs. (|3.2|) and Eq (3.8), the evalua- 
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formalism by defining an effective spectral density S'off(w) 
which should incorporate essential features of this broad- 
ening without simulating specific features. Having this in 
mind, an effective Gaussian density of states 5„(w) is as- 
sociated to each level \n) of energy En 



jhLO - Enf 

2(?ia„)2 



27r(T„ 



(3.14) 



The spectral effective density Scs{ujnm) associated to the 
transition m <— n is obtained by the convolution of iS„(aj) 
with Srn{uj) and with the excitation source spectral den- 
sity S{u!). For a monochromatic source the latter is 
also taken as a Gaussian centered at the modulation fre- 
quency that once integrated over all frequencies yields the 
square of the intensity of the modulation source. The net 
result is that Scs{iOnm) has the form 



5'off(w) = Soe 



(3.15) 



with Weff determined by the modulation frequency fl and 
the energies En and Em ■ The effective width aes contains 
information about the frequency widths of the excitation 
source and those of each level. 



IV. COMPARISON WITH EXPERIMENTAL 
RESULTS. 

We have tested the procedure described in last section 
to model parametric excitation in a specific experiment 
conducted at LENS. In this experiment "^^K fermionic 
atoms are trapped in a ID lattice, realized retroreflect- 
ing linearly polarized light obtained from a single-mode 
Ti:Sa laser at A=787nm, detuned on the red of both 
the Di and D2 transition of potassium, respectively at 
769.9 nm, and 766.7 nm. The laser radiation propagates 
along the vertical direction, to provide a strong confine- 
ment against gravity. The laser beam is weakly focused 
within a two-lens telescope to a waist size wq — 90 /zm, 
with a Rayleigh length zji=3 cm; the effective running 
power at the waist position is P=350mW. 

The trap is loaded from a magneto-optical trap 
(MOT), thanks to a compression procedure already de- 
scribed in IgI, with about 5 x 10^ atoms at a density 
around lO^^m""^. The typical vertical extension of the 
trapped atomic cloud, as detected with a CCD camera 
(see Fig. |]), is 500 /Ltm, corresponding to about 1200 oc- 
cupied lattice sites with an average of 400 atoms in each 
site. Since the axial extension of the atomic cloud is much 
smaller than zji, we can approximate the trap potential 
to 

_2r2 

V{r,z):^Voe ^0 cos^^kz); fc = 27r/A , (4.1) 



thus neglecting a 5% variation of Vq along the lattice. 
The atomic temperature in the lattice direction is mea- 
sured with a time-of-flight technique and it is about 
bOfiK. 

In order to parametrically excite the atoms we modu- 
late the intensity of the confining laser with a fast AOM 
for a time interval T~100ms, with a sine of amplitude 
e=3% and frequency The variation of the number 
of trapped atoms is measured by illuminating the atoms 
with the MOT beams and collecting the resulting fluo- 
rescence on a photomultiplier. In Fig.^ the fraction of 
atoms left in the trap after the parametric excitation is 
reported vs the modulation frequency f2/27r. Three res- 
onances in the trap losses are clearly seen at modulation 
frequencies 340 kHz, 670 kHz and 1280kHz. By identify- 
ing the first two resonances with the lattice vibrational 
frequency and its first harmonic, respectively, we get as 
first estimate luq ~ 2tt x 340 kHz. As we will show in the 
following, these resonance are actually on the red of loq 
and 2wo, respectively, and therefore a better estimate is 
Wo — 27r X 360 kHz. Therefore the effective trap depth is, 
from Eq. (9), Vq — 185 ~ 10.5 Tiwo. Since the atomic 
temperature is about Vo/3.5, we expect that most of the 
energy levels of the lattice have a nonnegligible popu- 
lation and therefore the anharmonicity of the potential 
could play an important role in the dynamics of paramet- 
ric excitation. Note that the third resonance at high fre- 
quency, close to 4wo, is not predicted from the harmonic 
theory. It is possible to observe also a much weaker reso- 
nance in the trap losses around 1.5 kHz, which we inter- 
pret to be twice the oscillation frequency in the loosely 
confined radial direction. Anyway, in the following we 
will focus our attention just on the axial resonances. 

As discussed in the previous Section, the overall width 
of the excitation assumed for our model system could 
play an important role in reproducing essential features 
of experimental data. Since the source used in the ex- 
periment has a negligible line width, it is necessary to 
model just the broadening of the atomic resonances. The 
spread of the transition energies due to the axial anhar- 
monicity is reported in Table I, while the broadening of 
each energy level, due to the periodic character of the 
sine potential, is estimated using Eq. (^). We now note 
that the ID motion assumed in Section II is not com- 
pletely valid in our case, since the atoms move radially 
along a Gaussian potential. Since the period of the radial 
motion is about 500 times longer than the axial period, 
the atoms see an effective axial frequency which varies 
with their radial position, resulting in a broadening of 
the transition frequency. Other sources of broadening 
are fluctuations of the laser intensity and pointing, and 
inhomogeneities along the lattice. We note that also elas- 
tic collisions within the trapped sample, which tend to 
keep a thermal distribution of the trap levels population, 
can contribute to an overall broadening of the loss res- 
onances. Since it is not easy to build a model which 
involve all these sources, we introduce semiemp irical ly 
an effective broadening for the r-th level (see Eq. (3.14)). 



4 



Recognizing that the width could be energy dependent 
we considered the simple expression 



(4.2) 



for several values of the constants Ai,Ao and p. When 
p — 0, i.e. for a constant value of the band width we 
were not able to reproduce the general experimental be- 
havior reported in Fig.|^. The best agreement between 
the simulation and the experimental observations is ob- 
tained for Ao = 0.0002, Ai = 0.0135, in units of uj^, 
and p ~ 3. Similar results are obtained also for slightly 
higher (lower) values of Ao,i together with slightly higher 
(lower) values of the power p. In Table I the resulting 
widths are shown for the lower twelve levels. Note that 
we have intentionally excluded levels 11, 12 and 13 from 
the calculation, since their intrinsic width is so large that 
the atoms can tunnel out of the trap along the lattice in 
much less than 100 ms Anyway, the inclusion of these 
levels proved not to change substantially the result of the 
simulation. 

The comparison of experimental and theoretical results 
is made in Fig. ||; the abscissa for the experimental data 
has been normalized by identifying uq with 27r x 360 kHz. 
As already anticipated, the principal resonance in trap 
losses appears at ~ 1.85 wq- This result follows from 
the fact that the excitation of the lowest trap levels is 
not resulting in a loss of atoms, as it would happen for a 
harmonic potential. On the contrary, the most energetic 
atoms, which have a vibrational frequency smaller than 
the harmonic one, are easily excited out of the trap. The 
asymmetry of the resonances, which has been observed 
also in j|] , is well reproduced in the calculations and it is a 
further evidence of the spread of the vibrational frequen- 
cies. The first interesting result obtained by our study 
of parametric excitation is therefore the correction nec- 
essary to extract the actual harmonic frequency from the 
loss spectrum. For the specific conditions of the present 
experiment, we find indeed that the principal resonance 
in the trap losses appears at Q ~ 1.85 t^o- Anyway, the 
calculation shows that the resonance is nearby this po- 
sition for all the explored values of Ao,i and p also for 
deeper traps, up to Vo=25huJo, and therefore it appears 
to be an invariant characteristic of the sinusoidal poten- 
tial. 



The result of the numerical integration of Eqs. ( 3.12 ) 
reported in Fig.^ reproduces relatively well the subhar- 
monic resonance, which in the harmonic case would be 
expected at loq. On the contrary, both experiment and 
calculation show that the actual position of the resonance 
is f2 ~ 0.9 luq. It must be mentioned that the accuracy of 
these result s is r estricted by the finite difference charac- 
ter of Eqs. ( 3.12| ) and by the fact that some noise sources 
which have not been included could be resonant at a 
nearby frequency. In particular, a possible modulation 
of the laser pointing associated to the intensity modula- 
tion is expected to be resonant at = ujq in the harmonic 



problem and it could play an analogous role in our 
sinusoidal lattice. 

The higher order resonance around 3.5 t^o observed in 
the experiment is also well reproduced by the calcula- 
tions based on first order perturbation theory. Note that 
a simpler approximation to the confining potential by a 
quartic potential Vq{z) = k2z'^ + k4z'^ would yield a reso- 
nance around 4wo and not 3.5 wq- Anyway, it is possible 
to understand qualitatively one of the features of this 
resonance considering a quartic perturbation of the form 
e{t)VQ to a harmonic potential. In this case the ratio of 
the transition rates at the 2ujq and 4a;o resonances is set 
by Eqs.O to 



|T(n±2,n)|V|r(?i±4,n)|2 cx 



(4.3) 



This result can qualitatively explain the absence of the 
corresponding high order resonance in the radial excita- 
tion spectrum (see Fig.^: since the radial trap frequency 
is a factor 500 smaller than the axial one, the relative 
strength of such radial anharmonic resonance is expected 
to be suppressed by a factor (500)^. In conclusion, high 
harmonics resonances, which certainly depend on the ac- 
tual shape of the anharmonic potential, are expected to 
appear only if the spring constant of the trap is large. 

In Fig.^ theoretical and experimental results for the 
evolution of the total population of trapped atoms at the 
resonant exciting frequency Q = 2wo are shown. Al- 
though there is a satisfactory agreement between the 
model and the experiment, we notice that experimental 
data exhibit a different rate for the loss of atoms be- 
fore and after 100 ms. This change is probably due to a 
variation of the collision rate as the number of trapped 
atoms is modified, which cannot be easily included in the 
model. The comparison of the experimental evolution of 
the trap population with and without modulation shows 
the effectiveness of the excitation process in emptying the 
trap on a short time-scale. 

We have also simulated the energy growth of the 
trapped atoms due to the parametric excitation, which 
is reported in Fig.|^. Our calculations show nonexponen- 
tial energy increase in contrast w ith what expected in 
the harmonic approximation, Eq. (3.7). The fast energy 
growth at short times is related to the depopulation of 
the lowest levels, which are resonant with the 2a;o para- 
metric source. The saturation effect observed for longer 
times is due to the fact that the resonance condition is 
not satisfied for the upper levels so that they do not de- 
populate easily. 



V. CONCLUSIONS. 

We have studied both theoretically and experimentally 
the time evolution of the population of atoms trapped in 
a ID sinusoidal optical lattice, following a parametric 
excitation of the lattice vibrational mode. In detail, we 
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have presented a theoretical model for the excitation in 
an anharmonic potential, which represents an extension 
of the previous harmonic models, and we have applied 
it to the actual sinusoidal potential used to trap cold 
potassium atoms. The simulation seems to reproduce 
relatively well the main features of both the spectrum 
of trap losses, including the appearance of resonances 
beyond 2a;o, and the time evolution of the total number 
of trapped atoms. 

By comparing the theoretical predictions and the ex- 
perimental observations the usefulness of a parametric 
excitation procedure to characterize the spring constant 
of the trap has been verified. Although the loss reso- 
nances are red-shifted and wider than what expected in 
the harmonic case, the lattice harmonic frequency can be 
easily extracted from the experimental spectra, to esti- 
mate useful quantities such as the trap depth and spring 
constant. 

We have also made emphasis on the need of modeling 
the broadening of bands with a negligible natural width 
in order to reproduce the observed loss spectrum. In a 
harmonic model this broadening is not necessary since 
the equidistant energy spectrum guarantees that a sin- 
gle transition energy characterizes the excitation process. 
We think that most of the broadening in our specific ex- 
periment is due to the fact that the actual trapping po- 
tential is not one-dimensional, and also to possible fluc- 
tuations and inhomogeneities of the lattice. 

To conclude, we note that the dynamical analysis we 
have made can be easily extended to lattices with larger 
dimensionality, and also to other potentials, such as 
Gaussian potentials, which are also commonly used for 
optical trapping. 
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5.944 
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7.513 
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8.240 


0.680 


0.082 


10 


8.920 


0.624 


0.092 


11 


9.544 


0.551 




12 


10.095 


0.402 




13 


10.497 







TABLE I. Energy spectrum in units of hcuo obta ined 
from the diagonalization of the Hamiltonian Eq. ( |2.4[ ) for 
Vo = 10.5 fiwo in a harmonic basis set with the lowest 40 
functions. The third column shows the band widths ar, 



Eqs. (3.14)and (4.2), used in the numerical simulations re- 
ported in Section IV. 
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FIG. 1. Absorption image of the atoms in the optical lat- 
tice, and shape of the optical potential in the two relevant 
directions. 
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FIG. 2. Experimental spectrum of the losses associated to 
parametric excitation of the trap vibrational modes. For the 
low and high frequency regions two different modulation am- 
plitudes of 20% and 3% respectively, were used. 



FIG. 4. Theoretical (continuous line) and experimen- 
tal (triangles) results for the evolution of the population of 
trapped atoms at the resonant exciting frequency 57 — 2ijQ. 
The circles show the evolution of the population in absence 
of modulation. 
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FIG. 3. Experimental (circles) and theoretical (lines) frac- 
tion of atoms left in the trap after parametric excitation vs 
the modulation frequency. The continuous line corresponds 
to the numerical integ ration of the first order perturbation 
theory equations (3.11) and the dashed line to the numerical 
integration of the finite difference second order perturbation 
theory equations (3. J). 
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FIG. 5. Calculated evolution of the average energy of the 
trapped atoms during parametric excitation at f2 = 2ljo. 
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